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DERIVATION OF THE EQUATIONS OF MOTION OF A SYMMETRICAL 
WING-TIP-COUPLED AIRPLANE CaNFIGUEATiaN WITH 
ROTATIONAL FREEDCM AT THE. JUNCTURES 
By Allaert A. Schy 

. SOMMARY 


The method of Lagrange multipliers Is used to take account of the 
dynamic effects of the ccmstralnts at the wing tips idien two Identical 
aliplanes are coupled to tHe' wing tips of a "mother" airplane. The 
resulting equations of motion of this symmetrical configuration are 
derived for one, two, or three degrees of rotational freedom at each 
Joint. 

The effects of aerodynaralc interference are ignored in the origi- 
nal derivation. It "is then shown hovr the' aerodynaiaic interference terms 
may he included when they are known. As an exajo^le, the interference 
terms arising from variations of the rolling velocities and angles of 
attack of the individual aiiplanes are treated. 


INTRODUCTION 

Wing- tip to wing-tip coupling of airplanes is being investigated t,s 
a Trtf>aTifl for. carrying fighter protection on bomber missions or for in- 
flict refueling. Because of the increase in effective aspect ratio in 
this configuration^ the outer ^parasite" airplanes can probably be 
carried more efficiently than in any other manner (references 1 to i^-). 
The method of tandem coupling does not have this advantage of increased 
aspect ratio and has, in addition, proven to be inefficient because the 
fear airplane must fly in the downwash of the front airplane. 

The structural loads on the wings in the wing-tip-cpupled configu- 
ration may be minimized by allowing rotational freedom at the joints. 

In this paper it is shown how the method of Lagrange multipliers may be 
used to analyze the dynamic effects of the constraints at the wing tips. 


2 


NACA EM L51G12 


The resulting equations of motion for small disturbances of a symmetrical 
configuration from its trim condition are derived for one^ two, or lihree 
degrees of rotational freedom at each joint. Because of the symmetry of 
the configuration, the equations may he separated into independent 
lateral and longitudinal modes. 

The primary purpose of this investigation is to analyze the purely 
dynamic effects of wing-tip coupling; therefore the equations of motion 
are first deidved without considering the aerodyneunic interference 
between the airplanes. Such Interference effects, hc5wever, may.be 
Important. Therefore, the types of aercjdynamlc interference terms 
which may occur are discussed from a general point of view, and it is 
shown how these terms, when they are known, can be included in the equa- 
tions of motion. As an example, it is shown how the Interference terms 
arising from variations of the rolling velocities and angles of attack 
of the Indivldiml airplanes may be Included in the equations of motion. 
These terms are believed to be the most Important Interference effects 
in this type of-coupllng. 

The effects of aeroelastlclty are ignored in this discussion. 


SYMBOLS 


X, Y, Z 



conventional stability axes for describing airplane 
motions; zom^onenta of aerodynamic forces along these 
axes; also con^onents of fixed positions of coupling 
joints in these axes 

steady-state velocity^ taken along the steady-state 
X-axis 


t 


time 


X, y, z con^onents of disturbance displacement of airplane along 

X-, Y-, and Z-axes, respectively 

u, V, V components of disturbance translational velocity along 

X-, Y-, and Z-axes, respectively (dx/dt^ dy/dt^ dz/dt 
taken in stability axes) 

u*, a nondimens ional forward velocity, sideslip angle, and 

angle of attack, respectively 

components of disturbance rotation about X-, Y-, and Z- 
axes, respectively 


0, 0, t 


KACA EM L51G12 


3 




an. arbitrary factor (Lagrange’s "undeteimlned multipliers”) 


5 


Indicates a small virtual inclement 


i index denoting degrees of freedom of a mechanical system 

q.^ generalized coordinates used to describe a mechanical 

system; also used^ in paxticular, to denote degrees of 
freedom of the three airplanes 


^i 

Ei 

m 

W 

^xz 

r 

L* 

h 


S 

P 



^Z 

L, N 


generalized applied force in airplane degree of freedom 

sum of the inertial, vei^t, and aerodynamic forces in 
unconstrained equations of motion of an airplane 

mass 

• veig^t 

moments of inertia of airplane about the X-, Y-, and Z- 
axes, respectively 

product of inertia of airplane 

flight-path angle 
lift force 
wing span 
wing chord 
wing area 
air density 


components of applied forces along stability axes 

conponents of applied moments about X-, Y-, and Z-axes, 
respectively; also cojtponents of aerodynamic moments 


angle of filter deflection about hinge, when hinge- type 
coupling is used 


a 



h 
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aiigle "bet-ween hinge axis and steady-state X-axis -when ’ 
hinge axis is parallel to steady-state XZ-plane 




angles between hinge. axis and the steady-state X-, Y-^ 
and Z-axes^ respectively, when the hinge is arbitrarily 
oriented 


^f^f 


SfCf 


S^c^'’ 

Sc3 


Js. A. 

Vq <it 


Vo dt 
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sq.uare of the nondlmensional radii of gyration in roll, 
pitch, and yav, respectively 

\fiLb^ mc^ mh^/ 

/-^xz\ 




^^b = 


^^c = 


^ “ 

Cx 

Cy 

Cl 

Cm 


nondimensional product-of- inertia parameter 


m 

pSb 

3U 

pSc 

^ <=0® ’'o 


longitudinal- force coefficient 
lateral- force coefficient 

\qSj 

normal-force coefficient 
rolling-moment coefficient 


(*) 


pitching-moment coefficient 1 

\qsc/ 

yawing- moment coefficient (-§=-l 

\q.Shy 


SCr 


= 


i^) 


Sc-, 


n-n “ 


'ft) 
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p • • • 

Whenever a, u', v, 3, w, v , a, p, 6, and t|t are used as suTd- 

SY 

scripts, a derivative is indicated. For example, Yy = !^ and Gy = r— ±-. 

ov P dp 

Kote: Unconventional stability derivatives caused by aerodynamic inter- 

ference are distinguished by parentheses; namely^ 



aerodynamic Z- force on "bomber due to unit filter angle of 
attack and unit flgjhter rolling velocity, respectively 


aerodynamic rolling moment on "bom'ber due to unit fighter 
angle of attack and unit filter rolling velocity, 
respectively 

{^t\> (Zf)^ 

aerodynamic Z-force on filter due to unit "bom'ber angle 
of attack and unit-bomber, rolling velocity, respectively 


aerodynamic rolling moment on the flgfhter due to unit 
bom'ber angle of attack and unit bomber rolling velocity, 
respectively 

5)f 

aerodynamic Z-force on filter due to unit rolling 
velocity of same filter 


mCA EM L51G12 


7 



aerodynamic rolling moment on filter due to unit angle 
of attack of same filter 


Su*bscripts: 


0 
f 

1 , 2 

1 

k 


Bf a 


initial conditions : 

filter parameters 

parameters of . left and rlg^t filters, respectively 

index denoting the degrees of freedom of a mechanical 
system 

index denoting a particular equ?.tlon of condition con- 
necting the degrees of freedom, or the corresponding 
undetermined multiplier 

eymmetrlcal and antisymmetrlcal components of comhined 
filter motion, respectively (see equations (9)^ (lO)^ 
and (l6)) 


Note: Ho subscript Is used on bomber parameters. 


ERELIMINAET DISCUSSIOH 


The configuiation to be considered consists of three alip)lanes 
coupled vlng tip to ving tip. The two outer airplanes are assumed to 
be exactly alike. Since the central airplane Is carrying the two outer 
airplanes along as parasites, these airplanes will probably be smaller 
than. the central airplane. For convenience, the central airplane will 
he referred to as a bomber and the outer aliplanes will be referred to 
as fighters, althou^ the analysis will actually be q.uite general and 
the only real restriction will be that the outer aiiplanes, including 
their coupling to the central airplane, be exactly alike. 

The Oldinary equations of motion for the three airplanes flying 
independently will be modified by the Interactions between the air- 
planes. These Interactions are of three types: the dynamic interactions 

arising from the coupling at the wing tips, the aerodynamic interactions 
arising from the interference effects on the air flow over the airplanes 
due to their proximity, and the elastic interactions. In the present 
case the airplanes are assumed to be rigidj therefore, elastic effects 
may be ignored. 
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The primary puip>ose of the present paper is to take account of the 
dynamic effects of the constraints at- the wing tips and set up a method 
for- obtaining the equations of motion of the entire configuration for 
small disturbances from a trim condition. This method will first be 
developed without considering the aerodynamic interference effects. 

The equations so derived would be correct for airplanes coupled together 
with a large enough gap between the wing tips so that aerodynamic inter- 
ference could be ignored. Later, it will be shown how the effects of 
aerodynamic interference may be considered without any essential modi- 
fication of the method of obtaining the equations of motion of the 
configuration. 


DERIVATION OF EQJJATIONS OF MOTION 
WITHOUT AERODYNAMIC INTERFERENCE 


In order to describe the motion of three independent airplanes, 
eighteen degrees of freedom consisting of the ordinary six degrees of 
freedom for each airplane must be considered. If the wing tips are 
assumed to be connected, then the tiranslational motion of adjacent wing 
tips must be the same. The translational motion of the wing tips at' 
each point of connection is expressed in terms of the degrees of free- 
dom of the individual aiiplanes. Then the expressions for the trans- 
lational velocities of adjacent wing tips may be set equal to each 
other. Thus, for each connection, three equations are obtained relating, 
the degrees of ,freedom of the aiiplanes, since one equation is obtained 
for each component' of the translational velocity of the wing tips. There 
are, therefpre, six equations of condition (or constraint) connecting 
the ei^teen degrees of freedom for the case of complete rotational 
freedom at each connection. The equations of motipn of the total system, 
taking account of the effects of the constraints, will be derived using 
Lagrange *s method of undetermined multipliers (references 5 and 6). 

The equations of motion for each independent airplane are referred 
to stability axes fixed in the- separate airplanes (see fig. l). For 
the present discussion it will be assumed that in the equilibrium con- 
dition of - the configuration the fighters may be trimmed, so that all six 
wing. tips lie on a stral^t line; then in the steady state all three 
sets of airplane axes may be taken parallel. For each airplane, the 
conponents of the wing- tip velocities along axes fixed in space parallel 
to the steady- state axes may now be written in te-rms of airplane degrees 
of freedom as follows: 


Vo + u + za - Yt 


(la) 
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V + Vot + - Z0 (lb) 

V - Vq© + - xe (Ic) 

• • • 

■where v, v and are the disturbance translational and 

rotational velocities of the airplane axes, Vq is the steady-state 
velocity along the X-axis, and X, Y, and ZL are the fixed positions 
of the -wing tip in airplane coordinates. 

The subscript f -will be used to indicate the filter parajaeters 
connnon to both filters. The filter to the left of the boaiber vill be 
indicated by using the subscript 1 on its variables and the filter to 
the ri^t of the bomber vill be Indicated by using the subscript 2, T^o 
subscript will be used for the bomber. The constraining conditions are 
that the velocity of the ri^t -wing tip of the left filter must equ4l 
the velocity of the left wing tip of the bomber end the velocity of the 
left -wing tip of the ri^t filter must equal the velocity of the ri^t 
■wing tip of the bomber. For exan^le, the t^wo equations giving the 
X-velocities of the points of connection are, from expression (la), 

Vo + u + Z0+|t=Vo + u^ + ^ 

Vo + a + Z0 - I t = Vo + + Z ^©2 + ^ tg 

Solving these equations and the similar equations of condition for the 
translational velocities of the filters gives 

u^ = u + . Z0 + ^ ^ 2 ^ ' ( 2 . 1 ) 

= V + V^ijr + xt - Z0 - - X^y^ + Z^l . (2. 2) ' 

= V - Vo0 - I 0 - X© + Vo©i - ^ 01 + Xf0i' ■ (2. 3) 

• "h * • • 

U 2 = u + Z© tJ’ - Z ^©2 - "^2. 



(2.4) 
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V2 = V + Voi + xi - Z0 - (2*5) 

V2 = w - Vo6 + I 0 - X0 + Vo02 + ^ ^2 ■*■ ^^2 f 2. 6) 


These equations are the equations of condition imposed by the constraints 
•yftien the axes of all the airplanes are parallel. Inasmuch as small angles 
occur between the sets of axes during the motion, equations (2) should be 
modified to be exact by introducing the proper trigonometric functions of 
these variable angles. However, it is apparent that equations (2) remain 
valid to first-order accuracy as long as the disturbances are small. 
Therefore, for the purposes of stability analysis these equations are 
valid equations of condition if the airplanes have been alined as previ- 
ously described in the trim condition. Actually, other possible trim 
conditions exist in which there may be steady-state -angles between the 
bomber and filter; therefore, the bomber and filter axes are inclined 
at an angle to each other. In such cases the proper trigonometric func- 
tions of the steady-state angles must be Introduced into the equations 
of condition. Since the presence of these constant factors does not 
fundamentally alter the method of obtaining the equations of motion, the 
discussion of such trim conditions will be deferred until later. 


Application of Lagrange *s Method of Undetermined Multipliers to 
the Case of Complete Rotational Freedom of the Wing Tips 

Lagrange's method of undetermined multipliers provides a convenient 
means of taking account of - the constraints in the present problem. This 
method is based on d'Alembert's Principle, ■vdalch is an extension of the 
principle of virtual work to dynamics. For the simple case of the motion 
of a set of n particles, d'Alembert's Principle states thg.t, for 
"virtual" displacements 8q^, 



(Qi - mj^q^)Sqj^ 


0 


where m^ is the mass of one of the particles, q^ is the displacement 
in the 1th degree of freedom, 'and is the corresponding applied 

force. In the present appllcaticai, let Indicate the sum of the 

inertial, aerodynamic, and welg^it terms in the unconstrained equation 
of motion corresponding to a given degree of freedom of an aliplane, 
le.t Fj_ indicate any additional applied force in that degree of freedom. 
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and let q.j^ indicate the displacement in that degree of freedom. Then 
d’Alembert’s Principle takes the form 


+ mrer - Y^v + mV^ir - i]fW sin cos 7o)&y ” 

+ IxzV + + . . . = 0 

( 3 ) 

Here 6y = 8q.y, 8+ s and so forth for all eigditeen degrees of 

freedom; and the factor of each Sq^ is the equation of motion in the 
corresponding degree of freedom, for the unconstrained airplanes, 

Lagrange’s method then requires the conditions relating the virtual 
displacements in the various degrees of freedom. These equations of 
condition can be obtained from equations ( 2) . For example, • if equa- 
tion (2.2) is multiplied by 6t, 


18 


i=l 


(T'i - Ei)dqi s 


8yy = 8y + 8t + X 8t - Z 80 - 8t - Xf 8t^ + Zy 80j^ 


Since the constraint is geometrical, 8t may be taken as zero (see 
reference 6, p. 58). With 8t = 0, each of the six equations of condi- 
tion obtained from equations (2) is multiplied by an- arbitrary param- 
eter (k = 1, 2, . . . 6) and the following equations result; 


Xl^ 8 x + Z 80 + I 5i}r - 8 xy - Zf 80i + ^ 81 ^ 1 ) = 0 (4.1) 


Xg^8y + X 8lf - Z 80 - 8y^^ - X^ 8tj_ + Z^ 80^^) = 0 (4.2) 


X^fsz - I 80 - X 80 - 8zy - ^ 802^^ + Xy 80y) = 0 


( 4 . 3 ) 
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+ Z 80 - I 
+ X 8 t - Z 

Xg^Bz + I 80 - X 


8 i - 8 x^ - Z^ 80g - = 0 (4.4) 
80 - 8yg - X^ 8i|f2 + Zf 50g) = 0 • (4.5) 
80 - SZg -f ^ 80g + X^ 802^ = 0 (4.6) 


Since these equations of constraint-are inte^rahle, 'six of the 
sixteen degrees of freedom can he eliminated hy using them (refer- 
ence 6 ) . It mi^t he noted that the exact equations of constraint, 
including the trigonometric functions of the variable deflection angles, 
■would not only lead to nonlinear equations of-motion hut- -would also lead 
to nonintegrahle constraints. In--fchis case none of the •varlahles could 
he eliminated and -tihe problem would he extremely complicated. 

Equations (4) are now added to equation (3), and the factors of 
each 8 qjL- are collected. The arbitrary may he chosen- such that 

six of -these factors should vanish, since there are six parameters X-j^. 

Since there are twelve independent variables in the system, the vari- 
ables associated with the remaining 'twelve factors may he considered to - 
he Independent and -the associated displacements 8 q^^ are therefore 

arbitrary. Since the displacements 8 q^^ are arbitrary, -the factors of 
the • remaining t-welve displacements 8 q^ must also vanish to satisfy the 

equation obtained hy adding equations (4) to equation (3). Thus the 
parameters may he chosen such that the factor of each of the 

el^teen displacements 8 q.j^ must -vanish and eighteen equations are 
obtained. For example, the factor of 8 y from equation (3) is 
(|-Fy + mv - Y^v + mVQTir - "tlW sin 7 ^ - 0 W cos 7 q), whereas from the sum of 

equations (4) the factor of . 8 y is X 2 ^ 5 * Set-fcing -the sum of these 
factors equal to zero yields the first equation. Kie ei^teen equations 
obtained in this fashion are the following; 


mv - YyV + mVQ- 4 r - ilrW sin 7 q - 0 W cos 7 q + X 2 + = Py 


(5.1) 
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-IVv + Iz’t - - -L^ ; + |(>.3^ - X4) + X(X,2 + X5) = N (5-2) 

-IyV - IxZ^ “ ■*■ ~ ^(^2 ■*■ ^5) ■*■ 2(^6 “ ^3) ” ^ (5.3) 

mu - X^u - X^jW + 9W cos 7 ^ + X-l + ( 5*^) 

-2^u + inir - - (Zq + mV^^Q + 9W sin 7 q + X^ + Xg = F2 (5-5) 

-My^u - M^v - M^w - M^e + lyV + Z(Xi + X4J - X^X^ + = M (5.-6) 

Vl - Tv/i + sin 7^ - 0iW^ cos 7^ - A-g = ^5. 7) 

+ ^ZfVi - - IxZf^l “ ^ x^^ - XfX2 = % (5.8) 

■'■ - ^^^1 ■*■ ^f ^2 " -^ ^3 = ^1 (5.9) 

SOS 7o - ( 5 . 10 ) 

- (Ze^ + Vo)®i + ’'o - ^3 = C5.11) 


-Mu^^l .- %^'Wi - - Me^0i + ZfX-1 + 


(5.12) 
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"^ v ^'^2 + " ^ Tjf ^2 " ^ XZ ^ i ^2 “ ^ 0^^2 " " ^ f ^5 “ ^2 


-' W ^ 2 . - ^ XZf ^2 " Hf .'^2 ^ ^ XJ ^2 “ + "f ^6 “ ^2 


“ f '^2 - ^^^2 - ^^^2 + ® 2 % ^ o ® 7 o - \ = 

X X • ^ 


“Zu^Ug + aifVg - Z^^2 - (2©^ + “f^o)®2 + ?'o “ ^6 = ^Zj 


- Mu ^'^2 “ ^^.^2 ■ *^^^2 " ^ 0+.®2 ■*■ ^ yJ ®2 “ ^ f^if ■*■ ^ f ^6 = ®*2 


(5.13) 
(5. lit) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 


In these equations the terms in the undetermined parameters represent 
the constraining forces arising from the wing-tip connections. Equa- 
tions (4) and (5) D-ow give twenty-four equations in the twenty-four 
unknowns consisting of the el^teen degrees of freedom and the six param- 
eters Xjj.. Moreover, since the constraints given by equations (4) are 

integrahle, the system can he con^letely described by only twelve inde- 
pendent degrees of freedom. The twelve Independent equations of motion 
can be obtained by first solving for the constraint parameters and 

then substituting these values into any twelve of equations (5). Any 
six of the original degrees of freedom may now be eliminated by using 
the constraint conditions in the convenient form given by equations (2). 

The most convenlent“method of carrying out this process is to 
obtain twelve linearly independent combinations of equations (5) which 
eliminate the parameters A^. This may be done in many equivalent ways. 
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*but th.e following twelve linearly independent combinations seem to be 
the simplest: 


Side- force equation: (5-1) 

Longitudinal- force equation: (5A) 

Normal- force equation: , (5-5) 

Bomber yaw equation: ^ (5*2) 

Bomber roll equation: (5*3) 

Bomber pitch equation: (5-6) 

Left- fighter yaw equation: (5«8) 

Left-filter roll equation: (5*9) 

Left- filter pitch equation: (5-12) 

Rig^t- filter yaw equation: (5-1^) 

Rl^t- filter roll equation: ( 5 - 15 ) 

Ri^t- filter pitch equation: (5-l8) 


- (5.7) + (5.13) 

(6.1) 

- (5.10) + (5.16) 

(6.2) 

- (5.11) + (5.17) 

(6.3) 

• X(5.1) + |[(5.10) - (5.16)] 

.(6.4) 

- Z(5.1) + ||C5.17) - (5.11)1 

(6.5) 

- X(5.5) - Z(5.1^) 

(6.6) 

(5.10) - 3Jf(5.7) 

(6.7) 

- Zj(5.7) - ^ (5.11) 

. (6.8) 

- Zf(5.10) + Xf(5.11) 

(6.9) 

- ^ (5.16) - Xf(5.13) 

(6.10) 

+ Zf(5.13) + ^ (5.17) 

(6.11) 

- Zf(5.l6) + Xf(5.17) 

(6.12) 


The six translational degrees of freedom of the filters may be elimi- 
nated from equations (6) by using equations (2). Then these are the 
twelve linearly independent differential equations in the twelve 
remaining independent degrees of freedom which determine the small 
motions of the configuration. Any other combination of equations (5) 
which eliminates the parameters may be written as a linear com- 

bination of these twelve equations. The first three equations give the 
combined translational forces on the system, the next three equations 
give the rotational moments on the bomber; and the final ’six equations 
give the rotational moments on the filters. The solution of this set 
of equations would give the translational motion of the bomber center of 


16 


KACA m L51G12 


gravity and the rotational motion of .the bomher and both filters. The 
translational motion of the filters is of little interest^ but may^ of 
course^ be obtained from equation (2), . . . 


SYMMETRIZATION OF EQUATIONS OF MOTION 


The equations presented are not in the proper form to show the 
symmetry of the system. The proper form can be obtained by expressing 
the fighter motions and forces in symmetrical and antisymmetrical com- 
ponents. For example, consider the rolling equations of the fighters. 
T]ie rolling equation of the left fighter, equation (6.8), is 

* ■V'Zf'i' + [^fVo - ^ - 





f * 

2 






01 - 0iZfWf cos 7o, - 


^(^o\ - ’'o)®i - (^xZf + 








{\ ^ ^f\Yo - Z#f ^ 


( 7 ) 


The rolling equation of the rig^t filter, equation (6.11), is the 
as equation (7) except that the subscript 2 replaces the subscript 
-b replaces b, and -b^ replaces b^j that is. 


same 

1 , 
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Vh^ 


^ - 


(^: 


/bb 


^ + ZfYv^)7ot + m^[-f - ZZf]0 + 


-’^s 


0 - 


bf. , b. 


^ V'" + “f ^ ^ T ® " ■# ■"' ^'^f " ^f)® ■" 


T Vo® + 




') h - Ht * ^fNf ■" ^f’^f)'" 


4-S 


02 - 02 ^fWf cos 7o + 02 ^ " 

^VqZv^ - ¥f sin 7o)®2 “ (^XZ^ + mfXfZfyi |?2 + "*^ 

- Z^Wf sin 7^t2 = I;2 + ^f^y 2 ^Z2 


^S' 


t2 + 


( 8 ) 


These eq.uations do not shofw* the symmetry of* the system^ since they con- 
tain both longitudinal and lateral degrees of freedom. These tvo equa- 
tions^ however, may be replaced by their sum arid their difference; 
physically, this substitution is equivalent to replacing the individual 
filter rolling moments by the equivalent symmetrical and antisymmetrical 
components of the combined fighter rolling momehts. The symmetrical and 
antisymmetrical filter rotations may be defied as follows: 


8 ^^ 1(61 + 62 ) ♦s =i(+i -♦ 2 ) (9) 

i»a ®a=|(®l-»2) ♦a - I - t2) (XO) 

Then, one-half the sum of equations (7) an<i (8) is 
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fZ/ - + %ZfV„ -.x(Lv^ + ^ zj* - 

(Zvj. + . 5^ z„^jl + 


' “f 


+ 


(¥ 


+ z. 




+ zJl + Z„Y \ + -f- Z 

0f Vf f Vf j 14 . V. 


- 


0aZfWf cos 7o - -f 0a + -#(^6^ + X^^.“ + -#(^o% “ 


.Wf ^o)®a - + Vf2f)^a ■" 4 




'^a 


[jL^^ + ZfY^^)Vo - ZfWf sin 7^ta = | [^^ + Lg + + F^^^ - 

t 

^fz, - %)] (11) 

Also, one-half the difference of equations (7) and (8) is 

Td 43 "bj:* • "bj? •• Tdj 5 / * fc 

-§ Zuf^ - “f / w + -f + mf X 0 + ^(nifVo + ZZ^^ - XZ^^)0 - 


-§■ %Vo0 + 




+ “f v~n“ 


0s - 


^f ■" ■" ^^^^f ) ■" 


i-"-f 


k - 0aZfWf cos 7o - V^f ^ e's + ^(Z0f + - Z^Z^^) 0 


^o)®s - (Ixz^ + “fXfZf^i + + Zf^vJ - 


4 ^f 


^s + [(Lvf + ZfXv^)Vo - ZfWf sin 7^+3 = ? - 


Lo + 


‘(^Yi - ^Yg) - + ^Zg) 


(12) 
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The variables occurring in eq.uation (ll) are v, ijr, 0, 00^, 

and ilTg^, while those occunring ‘in equation (12) are u, v, 9, 0g, 

6g, and i|fg. The set of variables occurring in equation (ll) may be 
designated as the lateral degrees of freedom of the entire system and 
those occurring in equation (12) may be designated as the longitudinal 
degrees of freedom. That is, the lateral degjrees of freedom of the 
system are the oiilnary lateral motions of the bomber and the anti- 
symmetrical congjonents of the combined filter motions as defined in 
equations (lO). The longitudinal degrees of freedom are the longi- 
tudinal bomber motions and the symmetrical coii 5 >onents of the filter 
motions as defined in equaticms (9). 

If the remaining two pairs of flg^hter equations (equations (6.7) 
and (6.10), and equations (6.9) and (6.12)) are similarly replaced by 
one-half their sums and differences, by maM ng the substitutions given 
in equations (9) (lO), the same separation of variables occurs. 

Finally, using equations (9) and (lO) in the six bomber equations com- 
pletes the separation of the twelve equations of the entire system into 
six lateral equations and six longitudinal equations, which are inde- 
pendent of each other. These equations are presented in appendix A. 

•The longitudinal and lateral equations may now be treated sepa- 
rately. The stability of the longitudinal motions may be detemnlned by 
the well-known method of e:s 5 >anding the determinant of the longitudinal 
equations and evaluating the roots of the resulting characteristic 
equation. The. same may of course be done independently for the lateral 
motion. Actual solutions may be obtained by any of the well-known 
methods for solving sets of linear differential equations, such as the 
Laplace transform method, some step-by-step method, or by use of an 
analog calculator. If the rotational motion of a given filter is 
desired, it is sln 5 >ly necessary to use equations (9) and (lO). For 
example, 0]_ = 0g + 0g^, and 02 “ 0a “ 0s* ^ older to obtain the indi- 
vidual filter motions, both the lateral and longitudinal equations 
must be solved, since these motions contain both the lateral and longi- 
tudinal modes of motion. 

These equations hold for the most general ts^pe of coupling likely 
to be encountered. In practical cases many sir?)llflcations will prob- 
ably be possible. For exsuirple, Z and will generally be much 

"h 

smaller than ^ and respectively. Also, for unswept wings X 

^ b.p‘ 

and will be much smaller than ^ and respectively. These 

facts in conjunction with the approximate magnitudes of the familiar 
stability derivatives of the individual airplanes show that it may be 
possible, in practical cases, to ignore many terms in these equaticms. 
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Note that the set of lateral equations is eleventh order in the 
time- derivative operator and the longitudinal equations are tenth 
order. Also^ neither characteristic equation has a zero root in the 
general case. The solution of these most general equations of motion 
may, therefore, he a formidable task. In order to determine the effect 
of varying parameters- on the motion of the configuration, it vould be 
necessary to use an analog computer. The filter rotational motions 
must be restricted in order to reduce the differential order of the 
problem. 


CONSIDERATION OF STEADY-STATE FIGHTER DEFLECTIONS 


In deriving the equations of motion it -was assumed that there vere 
no steady-state angles between the airplanes in the steady-state con- 
dition. ' Inasmuch as such angles may exist, their effect is now 
considered. 

In the practical case, it seems reasonable to assume that the 
yawing moment of the fighters due to their drag will be trimmed out so 
that the steady- state filter yaw angles vanish. Since the X-axis is 
arbitrarily chosen along the initial fli^t path, the Initial pitch 
angle also vanishes. The remaining angle which may be considered is a 
possible angle of steady-state "droop" of the filter wings. 

In the following analysis, 7 ^ is assumed to be zeroj then the 
steady- state roll angles are determined by the following equations of 
trim for side force, vertical force, and fighter roll about the left 
and right coupling joints, respectively: 


L'f sln( 0 i)^ + 3111(02)^ = 0 
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The first of these equations in 5 >oses the obvious condition that 

( 0 i)q “ “(^2)0' (^a)o “ (^2)0 “ 

and the last two equations are Identical. The remaining conditions then 
are 


li ^ cos 0 q 


L* = W + aff (1 - cos^^q) 


(13) 

Uk) 


These equations show that the configuration may he trimmed with the lift 
on each filter not sufficient to support its own wei^t, in which case 

the trim angle of droop will he = cos“"^ W^* Equation (l4)- also 

shows that the homher lift will have to he correspondingly Increased to 
carry the unbalanced wei^t. 

Note that the trim equations are satisfied for negative droop also, 
since 0^ appears only in the cosine function. This fact implies that 
there' is a trim condition with the' filter wings poised above the homher 
wing level.’ However, a sinple consideration of the effect of sll^t 
variations of filter lift or roll angle on the fighter rolling moment 
about the hinge shows that the upper trim position is statically 
unstable, the = 0 position is a position of neutral equilihlrum, 
and the lower trim position is statically stable. It might therefore 
be desirable to allow the filter wings to droop somewhat if the bomber 
wings are capable of sustaining the additional steady-state load. It 
should be noted that the aerodynamic interference effects have been 
l^ored in this discussion. These effects ml^t make positions of nega- 
tive droop stable. 

The case of steady- state filter roll angles, such as have just 
been described, may be handled by making the appropriate change in equa- 
tions (2.2), ( 2 . 3 )^ (2.5)\7 and (2.6). Equation (l) now inplies, when 

(1*4)0 - -(*'1)0 - 0 o, ttot • 

* * * * * I * 

V + + XMr - Z0 = (vi + - Zf03^)cos + 

(vi - ^ 01 - Xfei)sin 
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V + VoV 

+ - 

II 

("■2 ■ 

<■ Votg 4 

■ V2 - 

Zf^2)°°® ^0 - 




(Vg . 

■ ^o®2 “ 

^0.- 
2 >^2 

. X^e2)sin 0o 

w - Vo© 

-|fi 

- X© = 

(-1 

- Vi 

■bf> L 
1*1 

- X^e^)cos- 0O 




(^1 

+ Vo^t^i 

+ Vi - 

Zf0i)sin 0o 

w - Vo© 

+I0- 

- X© = 

(w2 

- ^0^2 

- -#02 

- X^Sg^cos 0o 


(^2 + ^ q '^2 ^2 " 

These equations give the modified equations of condition and replace 
equations (2.2), (2.3)> (2. 5)^ aJ3d (2.6) as follows; 

» -Vo+i - ^f^l + (-V + + xi - Z^jcos 0Q - 

(w - v^e - I 0 - Xe)sln 0O (15.1) 

Vg = -Votg - + (v + Vgi + xir _ z^Jcos jJq + 

(w - Vo© + I 0 - xejsin 0 (15.2) 

^1 = '^o®i - ^ ^1 + ^f®i + (v - Vo© - I 0 - Xe)cos 0o 


(v + Vot + XV - Z0)sin 00 


( 15 . 3 ) 
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Wg = 7^02 + ^ 02 + Xf©2 + (v - + I 0 - Xe)cos - 

(v + Vgt + xjr - Z^^sln (15A) 

It is easily seen that these conditions will not destroy the symmetry 
properties of the final equations of motion. Actually^ the most logical 
way to derive the equations of motion in symmetrical form would he to 
use the equations of condition in symmetrical form (sum and difference 
equations, with symmetrical and antisymmetrical filter variables), 
together with the Independent filter equations in symmetrical form 
replacing the usual independent filter equations given' in equations (5)* 
This method was not used only because it was not desired to confuse 
unnecessarily the discussion of the application of Lagrange* s method by 
introducing questions of symmetry. 

In symmetrical form, equations (15) become 

Xf+a + + (t + Voir + X* - Z0)cos + ^0 sin (l6.1) 

xA + Zf08 - - Yo® - xe)siJi 0O (16-2) 

■Ws = V-o0a - ^ 03 + XfSg + (v - Vo0 - Xe)co8 0 q 

Vg, = = Yo0a - + Xf0a - cos 0o + (■"■+ ^ - Z0)sin 0 q 

It can "be seen that Vg_ and ■which will enter Into the lateral 

equations of motion, ^contain only lateral variables, "whereas Wg and 
Vg, -which enter only into -fche longitudinal equations, contain only 
longitudinal variables. Therefore, the symmetry properties of -t^e final 
equations are preserved even when steady-state roll. angles are assumed. 
Physically it is clear that any steady-s-tate ang3.es -vdilch preserve the 
symmetry of the steady-state configuration -will give rise to separable 
lateral and longitudinal equations. 


(16.3) 

( 16 . 


'»'a s = - Va - 




^1 + ^2 ^ 
2 

vi - V2 


— ^ -^0+8 
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SIMPLIFIED EQUATIONS FOB MOEE EESTRICTED COUPLING 


Some of the possible sln5)llfication.a of the general equations of 
motion given In appendix A have already "been mentioned. For example. 

If any ' of the coupling- joint position coa^jonents" (X, Z, and Z^) 

should vanish, many terms In the general equations would not appear. A 
more fundamental simplification Is that In which the rotation at the 
coupling joints Is restricted to less than three degrees of freedom. 

In sudi cases additional equations of condition may he In^osedj two of 
the variables are thereby eliminated for each restriction In a ccm- 
pohent of the rotational motion. Thus, the number of equations and the 
differential order of the equations of motion which must finally be 
solved are reduced. 

‘ If the filters are restricted, by the design of the coupling, 'to 
rotate about only one or two axes Instead of being completely free to 
rotate, then the 'equations of condition will generally involve the 
direction cosines of the hinge axes ■with respect to the steady- state 
axes. The method of treating these cases may best be Illustrated by 
carrying out the case for rotation about one a,xls only, corresponding 
to an ordinary hinge coupling. 

The usual practical type of coupling Is a hinge type of connection, 
the axis of ■which is parallel to the symmetry plane of the bomber. If 
the hinge axis Is parallel to the XZ-plane, the fighters and bomber are 
rigidly connected in pf,tch, and two additional equations of condition 
may Immediately be -written as 


= 9 

(17.1) 

?= G 

- (17.2) 


Now if the X-axes In the Individual aliplsnes may be chosen such that 
the hinge axis- Is in the direction common, to the steady-state X-axes, 
then the additional two equations of restriction are simply 


(17.3) 

^^2 


( 17 . 4 ) 
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In general, such a choice of the steady-state X-axes would only he 
permissihle If body axes were used instead of stability ajces. ' In this 
case the, original equations of motion represented by in equation ( 3 ) 

would have to be written in terms of body axes. However, if stability 
axes are to be used in setting up the uncoupled equations of motion, as 
has been done in equations ( 5 )^ then the direction of the steady- state 
X-axes must be along the stea^-state velocity. If the hinge axis is 
assumed to be fixed in the bomber or in the filter wings, some steady- 
state angle will generally exist between the hinge .axis and the velocity 
direction. This angle will be called and is assumed positive when 

the hinge is pitched above the velocity in the stead.y state. Let the 
disturbance deflection of the fighter about the hinge axis be called a 
and be taken positive in the ' same direction as positive roll. Since the 
disturbance angles are small they may be treated as vectors, and the 
consonants of along the steady- state X- and Z-axes are 0^-0 and 

ilr^ - ijf, respectively. Therefore, 


01 > 0 = ai cos tIq 02 - 0 = (J2 cos tIq 

= -CT^ sin tIq tg - t = -CTg sin 


Eliminating and ffg from these equations yields 

(0]_ - 0)Bin tIq = (iff - f^)coe 
(02 - 0)sln tIq = - TlTg^cos 


Thus, the equations of condition corresponding to equations (lT-3) and 
( 17 . 4 ) are 


t + (0 - 0jtan Ti^ 
\lf2 = t + (0 - 0g)tan 


(17.3*) 
(17 A') 


In this case the hinge has been assumed to lie parallel to the 
XZ-plane. By using equations (17) the four variables 6^ and may 

be eliminated. In the most general case of hinge coupling, the hinge 
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axis will not necessarily Tae perpendicular to the Y-axis, hut will make . 
certain steady-state angles with all the axes. Let these angles he 
Ti^, and Ti^ with respect to the X-, Y-, and Z-axls, respectively. Then 
hy proceeding exactly as before, equations ( 17 ) hecome 



©I = e + 

©2 = 6 + ^{^2 
^1 *= t 

= lif + 


0) 


0 ) 


cos 

T 12 

cos 

’ll 

cos 

ri2 

cos 

’ll 

cos 

’I 3 

cos 

’ll 

COB 

’I 3 

cos 

’ll 


( 18 . 1 ) 


( 18 . 2 ) 


(18.3) 


(18.4) 


These are the general eq.uations of condition which must he used along 
with equations (2) in applying Lagrange’s method for hinge coupling when 
the hinge is skewed with respect to the symmetiy plane. 

In order to illustrate the modifications Introduced hy the addi- 
tional conditions, the simplest conditions, represented hy equa- 
tions (it. l) to (17.4), will he considered. By proceeding as before, 
the additional equations corresponding to equations (4) are 

Xrjfbe - 60^) = 0 (19.1) 

Xq(56 - S02) 

, 


= 0 ( 19 . 2 ) 

= 0 (19.3) 


0 


Xip(8t - 8+2) = 


(19. w 
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These conditions introduce the following additional terms into 
eq.uatlons (5): (^9 ^ 10 ) ^ equation ( 5 * 2 ), + X,g^ in equa- 
tion ( 5 * 6 ), in equation ( 5 . 8 ), -Xj in equation ( 5 . 12 ), “^10 

in equation (5.1^), and -Xg in equation (5.18). Six of the equations 
of motion given in equation ( 6 ) will not he changed and will therefore 
still he valid in this case. These are equations (6.1), (6.2), (6.3), 
( 8 . 5 ), ( 6 . 8 ), and (, 6 . 11 ) which are the translational equations and roll 
equations. The suhstitutions given hy equations (17) make it possible 
to eliminate four more variables, so that only eig^it variables remain. 
Therefore, when equations (17) are used in the six equations which 
remain unchanged, only two additional Independent equations are needed. 

In fact, only two more linearly independent combinations of equations (5) 
Tdiich eliminate the X^ are possible. In order to eliminate the addi- 
tional Introduced, the simplest combinations of equations would 

seem to be ~ 

(5.2) + (5.8) + (5.lii-) +^-^-^[(5.10) - (5.16)] (20.1) 

(5.6) + (5.12) + (5.18) (20.2) 


The equations of motion may now be symmetrized exactly as before. 
The nondimensibnal variables u*, 3, and a are used in place of u, 

V, and w^ respectively. Then, for the case of hinge coupling wibh the 
hinge axis along the steady-state X-axis, if^ it is also assumed that 
X = Z = Xf = Z^ =s 0 in order to obtain the sinplest possible case, the 

equations of lateral motion are 


(m + 2m^yrJ> - ^Yp + 2Yp^^p + (m + 2m^)vjr - ^(W + 2tff)sln 7 ^ - 
0 W cos 7o - 2 |^a^f ^ 

-(Kp + 2Np^)p + [iz 2Iz^ + ^f)^V - 1^ + 2t^^ + 

^(b + I- - 2Ixz^^ f 

- N + + 1*2 + 


( 21 . 1 ) 


( 21 . 2 ) 
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-LpP - Ixzt + 




2Vr 


t + (ix + “f .^)0 - 2^^'^ 


^ *0a - ^ Ztef0a = ^ • 


(21.3) 


"b.p (^"b + 'b-p') 

^ z,,.. - Isl 


W, 


^’f - H, 


, T^b^ V *bb.p , 


-Lp^P - + 

(%f + ”f - (iflf + ^ '^a.^K - i ^1 ^2 - %(\ - ^Zs) ■ 

(21.lt) 

and the equations of longitudinal motion are 

(ly + 2Iy^)e - (1% + 2M§^)I - + 2M^,^)u« - (Mi + - 

(Ma + 2M^)“ + ^ ^ “a/s = M + My + Mg (21.5) 

e(W + 2Wf)cos 7o + (“1 + 2mf)Voi* - (x^, + 2X^.^)u* - (X^ + 2X^^)a + 


^ ^^0S = ^X + ^Xy + ^Xg 


(21.6) 


-[z© + 2Zg^ + + 2my.)V^0 + 0(¥ + 2Wf)sln 7 q - (Zu* + + 

(m + 2m^)Voa - (z^ + 2Z^)a - ^ = F^ + F^^ + F^, 


2 

(21.7) 


^(Vo + Z 0 f)® “ 0 ^ sin 7o + ^ Z^t^u* - ^ V^a + ^ Z^^^a + 




^Xf + - ^0. •" W: 


'ayj^s - ij^l - ^2 - + ^Zg^ 


( 21 . 8 ) 
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It Ttii eht be noted that the equations of motion in this case could 
have been, obtained directly from the general case given in equations (6) 
by simply adding the three pitch equations and the three yav equations 
and using ijf = = ilf2 and 6 = ® ®2 obtain equations (20.1) 

and’ (20. 2). The reason for this is that the constraints in this case 
are the same as in the general case except that two additional con- 
straints are added. The method of introducing the additional was 

used in order to show how the nwre general constraints (such as given 
in equations (l8)) mig^t be handled. Also, if the general equations 
are not known, it is always more convenient to Introduce all the con- 
straints as done here, rather than first deriving the general equations. 

The system described by equations (21) is probably the sin^ilest 
wlng-tip-coupled ' system. The lateral equations can be seen to have a 
seventh-order characteristic equation, but, since there are two zero 
roots, the problem is essentially fifth order. The longitudinal motion 
has a sixth-order characteristic equation with one zero root and is 
therefore also essentially of fifth order. 

Clearly, there are many possible cases which are less general than 
that described by the general equations in appendix A but more general 
than that described by equations (21). It is hoped that the discussion 
of the modifications Introduced by additional constraints has made clear 
the methods for obtaining the equations of motion for these intermediate 
cases. 


NOKDIMENSICKALIZATION OF. EQUATIONS OP MOTION 


The equations of motion may be nondlmensionalized in several ways. 
The use of a method which brings in the conventional stability deriva- 
tives of. the separate airplanes would be desirable, however. One such 
method will be illustrated by applying it to equations (21). 


In order to nondlmensionalize the lateral "equations, equation (21.1) 
is divided by qS, equations (21.2) 6ind (21.3) a-re divided by qSb, and 
equation (21.4) is divided by qS^b^. Then, the following substitutions 


are made : 


Sfbf 
- “Sb“^ 




nondlmensional time-derivative operator, Ujj 
equations are 



JL A. 

Vo dt* 


and, for the 


The resulting 
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^2m^ .+ + 2kjp^^^Y + “ 

Tjr^C^^ + EAjC^^^tan 7 q - C^0 - 2A2C^^^a « Gy + A^^Cy^ ^Yg) 


(22.1) 


- 


-(% * ■" [Wz^’+ + Al(l * 

- A1A4C + I Aj(l + * 

I ^3(1 + ^k)%^ - I 0 Jl>bS< + >«^/xZj“li^a ■" 

I A3(1 + A^)C - A^Ai^C . 0^ + Ai(0„^ t <h^) * 

i * H){% - %) (22-2) 


-Oj^P + * ||a3(1 + ‘^k)<\,^ - CjJDb’f + 

I AiOz^B^fie. - Cj - I A 3 (Cz^- - ( 22 . 3 ) 



p + 2Alj%^^Kj[Zj.I>b^t 


J(l + A1j)C2^,^ - I AijOj D{,t + I A411^^D^20 . 


"■ “* 2 ) ■ i(\ ■ 


(22,4) 
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In order to nondimenslonalize the longitudinal equations, equa- 
tion ,( 21 . 5 ) is divided by qSc, equations (21.6) and (21.7) are divided 
by qS, and equation (21.8) is divided by qS^c^. Then, the following 

substitutions are made: Bn = Bo = pf^f ^ ^ ^ 

^ Sc * ^ gj.3 ' J S / c* 

•n *^f ”^f ■ 

and, for the nondimensional time-derivative operator, 
Kg * ^ The resulting equations are 

BlB5^n^^^ci2^s = Cj, + B^(C^ S) ^^^.5) 

+ 2BiHc^)DcU« - + 2B3Cx^,^^u» - 

+ ^B^CYo^Y + B3B5Cx^^Dc0s = Cx + B3(Cx;L %) 

^ ^I'^Zq^ + + 2830 ^^^) tan - 

+ 2B3 Cz^,^^u« + + 2 B 3 Cz^^a - 

+ ®3®5^2cx^°c0s “ C2 + B3^C2^ ■*■ ^Zg) ' (?2.?) 
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1 ^ ^ , 5516(1 + 

I =Za^)“ci «3 = i Bsfli - Ola) ■ ¥ + =22) 


COKSZDERATIOR OF AERODYNAMIC INTERFERENCE EFFECTS 


Becauae of the aerodynamic interference between adjacpnt wings ^ It 
is clear that the pressure distributions arising from given motions will 
be modified. Probably the most in^orteint interference effects are the 
changes in the lift distribution resulting from an angle of attack or, 
rolling velocity. From references 1 and 2 it can be seen that^ in addi- 
tion to modifications of the ordinary stability derivatives, certain 
unusual stability derivatives are introduced. These derivatives arise 
because the lift distribution lesulting from the motion of a given wing 
doeS‘ not vanish at the wing tip, but "spills" over onto the adjacent 
wing. This phenomenon gives rise to two ln^jortant dynamic effects. 
First, there are forces on a given airplane arising from motions of the 
adjacent one, and, second, there are coupling forces between the lateral 
and longitudinal modes of the filters. The second effect, arises from 
the asymmetry of the lift on the fighters. For example, a change in 
angle of attack on the fighter wing causes a larger change in the lift 
on the Inboard wing than on the outboard wing, because of the inter- 
ference effect of. the adjacent bomber wing. Therefore a rolling moment 
is produced. 

The purpose of this section is not to evaluate these Interference 
effects, but simply to show how they may be brou^t into the equations 
of motion when they have been determined by some theoretical or experi- 
mental method. In general, it' can be seen that the three types of 
Interference which have been mentioned in the discussion of angle-of- 
attack and rolling-velocity effects (namely, modifications of the ordi- 
nary stability derivatives, cross derivatives between the bomber and 
filters, and cross derivatives between the filter lateral and longi- 
tudinal modes) are, dynamically, the only types which can occur even in ' 
the most general disturbance motions. Therefore, the rolling- velocity 
and angle-of-attack interference effects will be considered in detail, 
and, if any other interference effects should turn out to be important, 
they may be treated in a similar manner. 
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For this discusslan. the modifications of the -mlues of the ordi- 
nary stahllity derivatives are unimportant, since they may he taken 
care of hy simply replacing the ordinary derivatives hy their modified 
values. The types of- unusual stahility derivatives will he distin- 
guished hy use of parentheses. For example, the homher rolling moment 
due to unit change of fighter angle of attack will he written (L)^^, 

the filter rolling moment due to unit homher angle of attack will he 
written and the filter rolling mcanent due to unit angle of 

attack of the same filter will he written 

The unusual stahility derivatives arising from filter angle of 
attack are Because of Ihe symmetry, the 

following definitions may he made: 


"(^)l " (^)2 “ (^)f 

(23) 

(Z)^ . (Z)„^ = (Z)„^ 

(2lf) 

- -(L)^ B (L)^ 

(25) 


Similarly, for the stahility derivatives associated with a homher 


angle of attack, 

i 

the. following definitions may he used: 



(^l)a - 


- 

.( 26 ) 


1 

H 

tl 


= (M„ 

(2?) 

For filter 

rolling-velocity derivatives 



11 


111 

( 28 ) 


-(Z)^^ = 

'">02 

III 

(29) 


- 

CM 

*«— s 


(30) 
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Finally, for 'bomber rolling- velocity derivatives, 



(31) 

(Li)^ - (La). = (Lf). 

(32) 


The forces associated with these derivatives may now 'be treated 
exactly as the ordinary aerodynamic forces were. These are simply . 
additional applied forces and may be considered in equation (3), in the 
proper equations of equations (5), or in the final combined equations 
of .motion. The Ijerms mentioned in equations (23) to (32) enter into 
equations (5-3), (5.5), (5.9), (5.H), (5.15), and (5.17). Now suppose 
that it is desired to find the additional terms in the final equations 
of motion arising from these unconventional stability derivatives. For 
siaipliclty, assixme the system is of the type described by equations (21). 
In this case the only final equations 'which are affected are equa- 
tions (21.3), (21.4), (21.7), and (21.8). 

In equation (21.3), the aerodynamic Interaction forces ■vriaibh must 
be considered may be added to the rl^t-hand side of the equation as 
additional applied forces. These additional terms are obtained with 
the use of equations (25), (26), (26), (30), and (31) and are as 
follows : 

- «2) + + fig) - |[(Zf)^= - (Zf)^a - i 


These tems have been written out in detail to make clehr the method of 
obtaining the additional terms. By the use of the equations of condi- 
tion for this case, ot^ and may be eliminated as follows: 



'2 “2 





= a ^ A _ ^ ^ 

2 4Vo ^ 2Vo 2 




~*^f- ^2 
2 Vq T 
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and 

4(“l * “ 2 ) - “ - ^ 5*= 

|(“1 - «2) = - 2^ f< - ^ ?>a (34) 

Tlie additional terms to "be ad.ded to tlie left-hand side of equation ( 21 . 3 )- 
are obtained by changing the signs of the previously derived additional 
applied forces and using equation ( 3 ^): 

- 4(2$)^^ + b(Zfy 


A similar process yields the correction terms for the other Idiree 
modified equations. When the proper collection of terms is made^ the 
additional terms from the rolling-velocity and angle- of- attack aero- 
dynamic interference forces on the system described by equations ' ( 21 ) 
are * 




b 0 





( 21 . 3 ') 



-2|(Z)aj -• (M^» ^ t 2 (Z )^^03 (21. 7') 


= In practice^ if the aerodynamic interference forces had been knovm, 
they vould have been included in equation ( 3 ) as aerodynamic forces in 
the associated -vdth each degree of freedom. Therefore^ they, would 

have appeared in the final equations in the proper , form. This more 
logical procedure was hot followed for two reasons: . First, no attempt 
has been made to determine all the interference forces which mi^t be 
important, since this paper is primarily concerned with the dynamic 
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rather then the aerodynamic effects of ving-tip coupling. Second, the 
consideration of aerodynamic interference in the midst of the e^qjla- 
nation of the dynamic effects would tend to confuse this explanation. 

The additional terms in equations (21. 3')^ (21.4*), (21.7*), and 
(21.8*) show that the aerodynamic Interference terms do not destroy the 
symmetry of the final eq.uatidns of the total system, since it will he 
noticed that the eidditions to the lateral equations (21.3) (21.4) 

contain o nly lateral degrees of freedom and the additions to the longi- 
tudinal equations (21.7) and (21.8) contain only longitudinal variables. 
The physical reason for this is as follows: Although the aerodynamic 

interference causes unsymmetrical pressure distributions on the fighters, 
yet with respect to the plane of symmetry of the total system (that is, 
the bomber plane of symmetiy) these pressure distributions are symmet- 
rical for symmetrical motions and anti symmetrical for antlsymmetrical 
motions. 


CONCLUDING KE34ARKS 


A method has been presented for deriving the equations of motion 
for the small disturbance motions of a symmetrical configuration of wing- 
tip- coupled aliplanes. Lagrange’s method of- undetermined multipliers 
was used to taJse account of the dynamic effects of the constraints at 
the -wing tips on the equations of motion for small disturbances of the 
unconstrained system of-aliplanes. 

By -this method, the equations of motion of a wing- tip- coupled 
system with three degrees of rotational freedom at each coupling joint 
were derived with the assunption of no steady-state trim angles between 
the fighters and the bomber and with the effects of aerodynamic inter- 
ference between the airplanes being neglected. It was shown that the 
twelve equations in twelve unknowns, which are necessary to describe 
this general system, may be separated into mutually independent lateral 
and longitudinal modes consisting of six equations eaqh. The lateral 
equations are eleven-th order in the time-derivative operator and the 
longitudinal equations are tenth order. The equations show the impor- 
tance of the position of the coupling joint in determining -the stability 
deri-vatives of the combined system. The slnple modifications caused by 
steady-state fighter deflections were also pointed out. 

The effects of restricting the ro-tatlonal motion at the coupling 
joints were then considered, and it was shown how the equations of motion 
could be derived for these cases with fewer degrees of freedom. The 
inportance of the effects of the orientation of the rotation axes on the 
equations of condition "in these more restricted cases was pointed out. 

The equations of motion were then derived for -the sinplest case of hinge 
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coupling with the hinge axis parallel to the steady- state X-axis and 
the position of the hinges lying on -the common steady-state Y-sxls of 
the airplanes. For this simplest case, the solution of either the 
lateral or longitudinal equations was shown to he a sixth- order • problem. 
A convenient method for* nondimensionalizing the equations of motion was 
illustrated hy applying it to these equations. 

Finally, it was shown that the aerodynamic interference forces 
could he treated in the same manner as the ordinary aerodynamic forces 
and would introduce no fundamental difficulties. 

Since the type of coupling with conplete rotational freedom would 
he most desirable from the point of view of decreasing the strains on 
the <5lng structures, it would seem advisable to investigate the sta- 
bility of this case by use of analog computers. In order to obtain 
accurate results the aerodynamic interference forces would first have 
to be evaluated by some theoretical or experimental means. Since the 
equations of motion are known, the computing machines could then be 
used to evaluate the effects of varying slgniflceint parameters, such 
as the positions of the coupling Joints, on the stability of the 
system. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPEKDIX A 

GENERAL EQUATIONS OF MOTION OF A SIMMETRICAL WING-TIP-COOPLED 
CONFIGURATION OF AIRPLANES WITH CCMPLETE ROTATIOML 
FREEDOM AT THE COUPLING JOINTS 

The lateral equations of motion are obtained as follows by making 
the conventional transformation to nondimens ional translational veloci- 
ties j that is, u’ s P s and a = where a and p are, 

*0 *0 'O 

to first order, the disturbance angle of attack and angle of sideslip, 
respectively: 

(m + 2mf)VoP - (Yp 2Yp^)p + + (mYo + ^ Yp^^Y - 


sin 7o 



«• 

(w cob 7o)0 + 2mfZf0a 

- 


- 2(Wf cos 

ro)9'a - + 

2X- 

Sf+a 



2{Ypj - 

sin 7o)^ 

= Fy + Fy^ + Fy^ 




(Al) 

s 

mXVoP + (XTp - 

Np)p + (12 

i + Y - (k^ - 

h mXYo 

+ X 1 ' 

ZTo ^ f> 

+ 


(xw sin 7 q)y 

- ixE? + (ii^ ^ + (™ 



L + 

bbf. ■bbf. „ • 

“f 2 2 Vq 


ZfXa^)ea + 



^bW^ cos 7 q - 

'%)«a = 

» - XFy + - 

S'Xg) 



(A2) 
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. iimj - (Lp + ZYp)p - + ^niOTo - ^ 2^,^^ - (ZW sin 7o)^ + 

(ix + Elf '*' “ 

^ Zu/a ^ 2a«/a ' + ^(^Of ’ ^fZ^’f ) -^ 

^Qk + T^(Za^ - Wf ro)0a = L + ZE^ + (A3) 


mfZfV^P - + ZfYp^)P + Z^,^ - XLp^ - XZf.Yp^) H 

“fZ^fV^t - + “f(~ip “ ZZf^jS + + 

- “ip Zlof')^ + ^ + ^Xf)*0a “ j^(^f%f 

-X)- 


+ -5 


?*a " (ZfWf.cos 7o)0a " (^XZ^ + 


+ _|- 2^, I - 


•- 


ta + 


“i^fZ-f) ^a 

(^f + ZlfYp^ - ZfWf sin ro)iIfe. - -f Xf^a^ + j^_^XfZ^ 
^fZ^’ ■*■ '^ 9 ^ ®s "' ^o^®a “ 

(I^ + Lg) ■+ Zf^Fy^ "^ ‘^^2) " '^(^^1 “ ^Z.2) 


1 

2 


(Al^) 
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-"iXfVoB * ”l(^ - 5=f)+ * - 

^ Xu'j.) - “AyJ* + - “Bf)* * ^(™Bf - 




Xa. 


4 


NA 

0f 




v^V“ir ^ " 


0a + (%^f ^°® ^o)0a + (^f + “f “f“)^*a '’' " 

^ "^a + (^Pf “ ^o)^a " 

b^ •• b 


“fZf -§ e*a + ^(ZfXutf - Xf^)Sa + -f(^f '^os Tq- “ X^^)9a = 


(% + ^2) + ^(^Xi “ ^Xg) " +■ ^2) 


(A5) 


^f^)0a “ ^a + “ ^'f)^a (^Xf ■*■ 

“f^f^ + DifZfS _ ^ ■*■ ” 

-’XfM^^ -H ZfM^,^) - XfZ^^ - ■ 

Hxj - ZfWf, cos >Q + XjWf sin 7o)®a = 2 " “2) " ^f(^l " ^^2) "*' 

“ %)\ 
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The longitudinal equations for this general case are 
(m + 2mf)VoU* - (x^t + 2X^,.^)u‘ - + 2X^)a + 2sr^' + " 

+ (W cos + SX^^Y + + “fMs - ^ JTb " 

2“M * 2(Wf ==■= r, - - 




(A7) 


.(Z^,, + 2Z^.^)u* + (m + 2m^)Yji - (Z^ + 2Z^)a - 2m^0 + [^(XZ^^ - 

ZZu*f) - “^'^o “ 2mfVo - Z^e + sin 7q + ^a^Y - “ft>f?s + 

< 

• r 

^ Za^s - ^ Zu«/a + ^ ^ f ^6 


®s + 


sin 7 q - ^a^Ye “ ^Z + ^Z^ + ^Z, 


(A8) 


-mZVoU« + (ZX^t - XZ,,, - M^,)u* + (mJCVo - M^,)i + (ZX^ - XZ^ - M^)a + 


1^0 - (Mg + XZg + mXVQ^e + (Xtf stn - ZW cos' 7 q)0 = 


(A9) 


M + XF 2 - ZPx 


HICVI 
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^ Z^,^u« - Of ^ Voi + Z^a + ^ Xe + j|^^ZZ^,^ - 


^■P 

“f •# V, 


^ ^ s) ^ 


z^e + liv + 




■*■ “f ItJ^s “ 


0g - (ZfWf COS 7o)0s - (^XZf + + 


+ ZfXf?p^ + 5^ + ZfTp^ - 


ZfWf sin 7 q 
b, 


l®s + 


)ts - “fXf -f e’s + j^(XfZ^3^ - ZfZut^) + ^ ^0^®! 

“|•(^ “ % ^o)®s 2 [^ " ^2) + (%2. " ^^2) “ 

^(^Zi ^ (AlO) 


“f ^ V - ^ ^ + “f I Z® - Z3Cu.^)0 


2 


■ (^% ■'■■”f?A')^s 




(^ V* Wsf - - 


^6 + (Vf =°= ’’o)i^B + ('Zj + “f*f^ + ”f + . 

^ ^'j.) - -ii + (%f - + ^f»f =1“ ro)ts - 

”f ^ Zf«B + - ^Aj,)«b ♦ ^(«f “« /o - ^)9s - 

Cl - Nz) ^ ^ %) - ^{\ - %) 


(All) 
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-iDfZfVoU* + - M^i^^u-* + ^nifXfVo - M^^^ct, + - 

- V^Zf - JnfXX^ye + j^^ZZfX^,^ - Z^XX^^ - 
XfZZ^«^ + XfXZ^^ - ZMut^ + - 




Id ^ ( 

“fZf ^ ts + ^(ZfXu.^ - XfZ^.^ - Mutf)ts + fef + “fZ/ + “fXf2 - 

^ %)®8 + |^M0f - % - + ^(XfZfXaf -*- 3CfZfZ^,^ - - 

zA*f + Z^.^ - + (ZfX^ - XfZ^ - ZfWf cos-7q + 

Vff ’’o - «aj.)®a - i[(“l ”2) - ^f(^Xj. ^ia) + *ffzi ♦ 


(A12) 
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Figure 1,- Stability axes. The motion of each individTial airplane is 
given in teirms of stability axes in that airplane. Arrows indicate 
positive direction, ’ ■ . 
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